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ABSTRACT 

We report the results of a long-term spectroscopic monitoring of the FS CMa type object MWC 728. 
We found that it is a binary system with a B5 Ve (Tes = 14000±1000 K) primary and a G8 III type 
(Tefi ~ 5000 K) secondary. Absorption line positions of the secondary vary with a semi-amplitude 
of ~20 kms“^ and a period of 27.5 days. The system’s mass function is 2.3x10“^ Mq, and its 
orbital plane is ~13-15° tilted from the plane of the sky. The primary’s usini ~110 kms“^ combined 
with this tilt implies that it rotates at a nearly breakup velocity. We detected strong variations 
of the Balmer and He i emission-line profiles on timescales from days to years. This points to a 
variable stellar wind of the primary in addition to the presence of a circum-primary gaseous disk. The 
strength of the absorption-line spectrum along with the optical and near-IR continuum suggest that 
the primary contributes ^60% of the U-band flux, the disk contributes '^30%, and the secondary 
^10%. The system parameters, along with the interstellar extinction, suggest a distance of kpc, 
that the secondary does not fill its Roche lobe, and that the companions’ mass ratio is q ~0.5. Overall, 
the observed spectral variability and the presence of a strong IR-excess are in agreement with a model 
of a close binary system that has undergone a non-conservative mass-transfer. 

Subject headings: stars: emission-line. Be; (stars:) binaries: spectroscopic; stars: individual: MWC 728 


This paper is partly based on observations obtained at the 
Canada-France-Hawaii Telescope (CFHT) which is operated by 
the National Research Council of Canada, the Institut National 
des Sciences de FUnivers of the Centre National de la Recherche 
Scientifique de France, and the University of Hawaii as well as on 
observations obtained at the 2.7 m Harlan J. Smith telescope of 
the McDonald Observatory (Texas, USA), 2.1 m of the Observa- 


torio Astronomico Nacional San Pedro Martir (Baja California, 
Mexico), 2 m telescope of the Ondfejov Observatory, Czech Re¬ 
public, and 0.81m telescope of the Three College Observatory, 
North Carolina, USA. 
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Fig. 1.— A 10^ X 10^ field around MWC 728 (marked by the 
circle). The comparison star (HD 281193) and check star (NO¬ 
MAD 1197—46760 labeled with a “C”) are shown. North is at the 
top and East is at the left. 

1. INTRODUCTION 

This paper continues a series devoted to studying ob¬ 
jects from the Galactic FS CMa type group consisting 
of ~70 members and candidates characterized by the 
B [e] phenomenon. Th e phenomenon was discovered by 
lAllen fc Swingsl (|1976ll and refers to the presence of for¬ 
bidden emission lines of [Fe ii] and/or [O i] in the spectra 
of B~type stars. In the vast majority of cases one also de¬ 
tects a strong exces s of IR radiati o n tha t is attributed to 
circumstellar dust. iLamers et al.l (1199811 recognized four 
subgroups of B[e] objects with known evolutionary sta¬ 
tus: pre-main-sequence Herbig Ae/Be stars (“HAeB[e]”), 
symbiotic binaries (“SymB[e]”), compact Planetary Neb¬ 
ulae (“cPNB[e]”), and supergiants (“sgB[e]”). They con¬ 
cluded that the B[e] phenomenon is found in objects at 
very different evolutionary stages but with similar con¬ 
ditions in their circumstellar envelopes. However, they 
were unable to c lassify ~50% of the objects selected by 
lAllen fc Swingsl (|1976l l and called them unclassified ob¬ 
jects with the B[e] phenomenon (“unclBje]”). 

The unclBfel object s were critically reconsidered by 
iMiroshnichenkol (|2007ll . Several possible evolutionary 
states (e.g., pre-main-sequence stars or supergiants) were 
rejected, and the group was renamed FSCMa objects. 
The main observational properties of this group include 
the following: 1) early-B to early-A type optical con¬ 
tinuum with strong emission lines of hydrogen, Fe ii, [O 
i], and sometimes of [Fe ii] and [O ill] (some absorption 
lines from the hot star atmosphere may be present as 
well, but they are frequently veiled by the circumstellar 
continuum); 2) a large IR excess that peaks at 10-30 
/rm and sharply decreases longward; 3) location outside 
of star-forming regions; and 4) a secondary companion 
(so far discovered in about one third of the group mem¬ 
bers) which can be a fainter and cooler normal star or a 
d egenerate object. 

iMiroshnichenko et al.l (120071) expanded the group with 



Fig. 2.— The spectral energy distribution of MWC 728. Loga- 
rithm of the flux normalized to that in the U-band is plotted vs. 
logarithm of the wavelength in microns. Symbols show the fol¬ 
lowing data sets: circles - optical and near- IR photometry from 
TableJ^ upward filled triangles - WISE data fl Wright et al.li20igh 
upwa rd open triangles - IRA S data, and filled squares - AKARI 
data JMurakami et ahl 120071) . The fluxes were dereddened using 
the interstellar extensi on law f r omISavage &: Mathi^ l|1979lh The 
dashed lines represent IKuruczl 1)199411 model atmospheres for the 
hot star with Teff=14000 K and the cool star with Teff=5000 K. 
The solid line shows total fluxes with a 60% contribution from the 
hot companion, 10% from the cool companion, and 30% from the 
circum-primary gaseous disk to the V-band flux. 

10 newly discovered objects found in the IRAS database 
by cross-identif ication with catalogs of op tical posi¬ 
tions. Later IMiroshnichenko et al.l (|2011f) reported 
^20 more candidates found in the NOMAD catalog 
()Zacharias et al.l l2005h using optical and near-IR color 
criteria. Although the origin of the B[e] phenomenon is 
understood at least qualitatively in all the groups with 
known evolutionary status, it still remains controver¬ 
sial whether FS CMa objects are close binaries recently 
evolved off main-sequence or a distinct group of proto- 


TABLE 1 

Summary of the spectroscopic observations 


Observatory 

Dates 

Range 

Resolution 

No. 

CFHT 

2004-2010 

3600-10500 

65000 

3 

OAN-SPM 

2005-2014 

3600-8000 

18000 

36 

McDonald 

2005-2013 

3600-10500 

60000 

11 

TCO 

2013-2014 

4250-7850 

10000 

56 

Ondfejov 

2012-2014 

6250-6750 

12500 

10 


Column information: (1) — Observatory name, (2) - range 
of observing dates, (3) - spectral range in A, (4) - spectral 
resolving power, and (5) — numbers of spectra obtained. 

Spectrographs used: ESPaDOnS lIDonati et al.l 1199711 at 
CFHT, REOSC at Observatorio Astronomico N acional San 
Pedro Martir (OAN-SPM), cs23 IITull et al.ll99dl at the Mc¬ 
Donald Observatory, Eshel from Shelyak Instruments at the 
Three College Observatory (TCO), and a Coude spectro¬ 
graph IlSlechta &: ^kodall2002ll at Ondfejov. 
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planetary nebulae Icf. lMirosliniclieril^l2007f) . 

The subject of the present paper is a relatively bright 
{V ^ 10 ma g) emission-line star MW C 728 that was first 
reported bv iMerrill fc Burwelll (jl949H . Some other com¬ 
mon identifications of the object are IRAS 03421-1-2935, 
HD 281192, and BD-|-29°611. A brief survey of obser- 
vations of MWC 728 found i n the literature is given by 
iMiroshnichenko et alJ (|2007D . who described its optical 
and near-IR spectrum and included it in the FS CMa 
objects group. They reported the following spectral fea¬ 
tures: double-peaked Balmer and [O i] lines in emission. 
He I 4471 A and Mg ii 4481 A absorption line from a 
mid-B spectral type star, and absorption lines of Fe i 
5328 A, Li I 6708 A, and Ca i 6717 A lines which indi¬ 
cated the presence of a cool star. We have been moni¬ 
toring MWC 728 for about 10 years using high-resolution 
optical spectroscopy and multicolor photometry to study 
its nature and evolutionary state. 

Details of the observations are provided in Sect. O 
Analysis of the optical spectra and behavior of the ab¬ 
sorption and emission lines is presented in Sect. [31 Sug¬ 
gestions about the system’s nature are discussed in Sect, 
m and conclusions are summarized in Sect. |5l 

2. OBSERVATIONS 

Optical spectroscopic observations of MWC 728 were 
obtained between December 2004 and February 2015. 
A summary of the facilities and instruments used in 
this project is given in Table |T] We have obtained 
116 spectra with a range of spectral resolving powers of 
R = 10000 — 65000 and signal-to-noise ratios from ~30 to 
>200. The spectra obtained at the McDonald Observa¬ 
tory (Texas, USA), Observatorio Astronomico Nacional 
San Pedro Martir (OAN SPM, Baja California, Mexico), 
Ondfejov Observatory (Czech Republic), and Three Col¬ 
lege Observatory (TCO, North Carolina, USA) were re¬ 
duced in a standard way with the echelle/slit package 
in IRAE0. Observations obtained at the C FHT were re¬ 
duced with the Upena and Libre-ESpRIT (iDonati et al.l 
[19^ software packages. Typical uncertainties of the 
wavelength calibration are <lkms“^ for the CFHT and 
McDonald data and ^l“2kms“^ for the other observa¬ 
tories data. Individual radial velocity measurements may 
have larger errors depending on the line strength and the 
profile shape. 

Photometric data for MWC 728 were obtained at the 
following sites. Three UBVRj observations were ob¬ 
tained at a 0.6 m telescope of the Crimean Astrophysical 
Observatory (CrAO) with a single-element photometer 
in 2008 and 2009 (f or the data reduction technique see 
iCrankin et al.l 1^00^ . BVRq observations of a 10' x 10' 
field around the object were obtained on 23 nights at 
aim telescope of the Tien-Shan Observatory (TSAO) 
of the Fesenkov Astrophysical Institute of the National 
Academy of Sciences of Kazakhstan between September 
2014 and February 2015. A 3056x3056 Apogee F9000 
D9 CCD camera with 12/im pixels and a 3-step cooling 
was used with a set of BVR filters. The TSAO data were 
reduced using Maxim-DL, while brightness of the object 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under a cooperative agreement 
with the National Science Foundation. 


TABLE 2 

Photometric observations of MWC 728 


JD 

V 

U-B 

B-V 

V-R 

J 

H 

K 

2450000-1- 








4748.3885 

9.79 

0.16 

0.32 

0.31 




4750.4548 

4753.5535 

9.75 

0.04 

0.23 

0.31 

8.28 

7.46 

6.90 

4755.5736 





8.49 

7.62 

6.74 

5140.9304 





8.32 

7.45 

6.76 

5145.5718 

5487.9896 

9.80 

-0.17 

0.25 

0.40 

8.32 

7.45 

6.66 

5488.9299 





8.38 

7.51 

6.73 


The UBVRj observations were obtained at CrAO and the 
JHK observations were obtained at OAN SPM. 


JD 2450000-1- 

V 

to 

1 

V-R 

6917.265 

9.80 

0.31 

0.22 

6918.269 

9.85 

0.29 

0.24 

6945.429 

9.76 

0.27 

0.23 

6948.458 

9.78 

0.28 

0.24 

6960.396 

9.85 

0.31 

0.26 

6963.485 

9.78 

0.31 

0.27 

6972.184 

9.78 

0.27 

0.23 

6984.312 

9.79 

0.27 

0.24 

6988.188 

9.79 

0.26 

0.24 

6993.251 

9.77 

0.28 

0.23 

7000.182 

9.75 

0.28 

0.22 

7001.170 

9.77 

0.28 

0.22 

7002.262 

9.75 

0.27 

0.22 

7003.297 

9.77 

0.27 

0.21 

7005.326 

9.80 

0.25 

0.24 

7006.295 

9.75 

0.27 

0.22 

7019.167 

9.78 

0.27 

0.22 

7020.149 

9.77 

0.29 

0.23 

7022.133 

9.81 

0.28 

0.24 

7023.242 

9.79 

0.28 

0.24 

7050.266 

9.75 

0.28 

0.22 

7058.242 

9.89 

0.33 

0.24 

7059.054 

10.01 

0.37 

0.28 


The BVRc observations obtained at TShAO in the 
Johnson-Cousins photometric system. 


and field stars were measured using the IRAF task imex- 
amine. The TSAO instrumental magnitudes were con¬ 
verted into the standard Johnson-Cousins BVR system 
using data for several open clusters. The TSAO pho¬ 
tometry was calibrated by comparing the brightness of 
MWC 728 to HD 281193 as a comparison star and NO¬ 
MAD 1197—46760 as a check star (see Figurcjl]). The 
brightness difference between the comparison and check 
star was found to be < 0.01 mag in all our images. 
Photometry of HD 281193 was obtained at CrAO and 
resulted in the following magnitudes: V = 10.00 mag, 
B — V = 0.28 mag, U — i?j = 0.18 mag. 

Near-IR JHK observations were obtained at the 
0.84 m telescope of the OAN SPM on 5 nights in 
2008-2010 with a 256x25 6 pixels IR camera Camila 
l| Cruz-Gonzalez et al.lll99^ . Th e images were calib rated 
by observing standard stars from iHunt et al.l 1)19981) . The 
near-IR brightness of MWC 728 in our data is ^20% 
higher than that in the 2MASS catalog recorded in De- 
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cember 1999 fcf. iMiroshnichenko et al]l2007f) . The re¬ 
sults for all the described photometric observations are 
presented in Table [2] The accuracy of all the data is 
0.02-0.03 mag. 

3. DATA ANALYSIS 

The spectral energy distribution of MWC 728 is shown 
in Figure [2] It is characterized by a slightly reddened op¬ 
tical continuum of a hot star and a strong IR excess pro¬ 
duced by circumstellar gas and dust. A detailed analysis 
of the optical absorption and emission lines is presented 
in the following subsections, while the components of the 
total flux are discussed in Sect.lH The dust component 
of the system will be analyzed in a follow-up paper. 

3.1. Absorption lines 

Absorption lines in the spectrum of MWC 728 can be 
divided into three groups. 

The first group contains relatively broad (FWHM 
^220 kms“^) and weak (I/Icont < 0.9) lines that man¬ 
ifest a moderate projected rotation of a hot star (see 
Figure[3l left panel). It includes He i lines (4009, 4026, 
4471, 5876, 6678 A), the Mg ii 4482 A, Si ii 6347 and 
6371 A, and the O i 7772-7775 A triplet. 

The Mg II 4482 A line has an average equivalent width 
(EW) of 0.22±0.03 A which is typically 1.7±0.3 times 
weaker than the EW of the He i 4471 A line. The EW ra¬ 
tio of these lines serves as an effective temperature (Teff) 
criterion. We used spectra of ^^60 B-type dwarfs and gi¬ 
ants as well as Be stars obtained at the Observatoire de 
Haute-Provence (OHP, R ^48000. [Moultaka et al.l[200li 
and TCO to calibrate the ratio (Figure[3l right panel). 
Since the EW of both lines can be affected by a contri¬ 
bution from the cool star and the circumstellar material, 
we need to explore this effect. 

As shown below in this Section, the cool star makes a 
negligible contribution in this spectral region. Any addi¬ 
tional circumstellar continuum would lower EWs of both 
lines. Therefore, removal of such a continuum would lead 
to a lower temperature estimate based on the Mg ii EW 
but to a higher one based on the He i EW. Although the 
lines may also be affected by radiation transfer in the 
disk, the net effect is probably small judging from the 
positions of Be stars in the right panel of Figure[3l The 
He I lines in the blue part of the spectrum of MWC 728 
(EigureS]) look weaker than those of BS 896 (Tes = 14000 
K) and comparable to those of Vega (Teg = 9500 K). 
However, the circumstellar contribution to the observed 
optical continuum estimated in Sect. 3] partially accounts 
for the He i lines weakening. With this analysis, we esti¬ 
mate the hot companion Teg = 14000T1000 K which cor - 
responds to a spectral type B5 (jCrav fc Corballvl 1 199411 . 
Additionally, comparison with spectra of B-type stars, 
which show no evidence for a noticeable amount of cir¬ 
cumstellar gas, suggests a projected rotational velocity 
of usinz ~110 kms“^ (see Figure[3l left panel). 

Atmospheric absorption wings are only occasionally 
seen in Ha, while higher members of the Balmer series 
do not fill their atmospheric profiles with emission. The 
latter exhibit deep and narrow line cores during peri¬ 
ods of the enhanced emission (starting from HJ, see Fi- 
gurc|4]). The Si ii lines are weak that is typical of dwarfs 
and giants, thus excluding a supergiant luminosity. The 


Paschen series lines and Fe ii lines are seen neither in ab¬ 
sorption nor in emission. They are most likely filled with 
emission components as well as veiled by the continuum 
of the cool companion and circum-primary disk. 

Most absorption lines of the hot star are variable, espe¬ 
cially He I lines which sometimes show emission compo¬ 
nents (see FigureE]). These components are seen in the 
line wings and probably due to the fast rotating material 
in the inner parts of the circum-primary disk. The detec¬ 
tion of only weak emission supports our determination of 
the primary’s Teg, because He i lines transition into full 
emission typically occurs at Teg ~ 18000 — 20000 K. 

The second group contains narrow lines of neutral met¬ 
als, such as Fe i, Ca i, O i, Ti i, and Li i (see Figurc[6]). 
They become clearly detectable at A ^5300 A and can 
be traced all the way to the red end of our optical spec¬ 
tra (^10000 A). The weakness of these lines is caused by 
a stronger contribution of the continuum from the hot 
star as shown below. It is not easy to constrain the tem¬ 
perature of the cool companion in this situation. The 
obvious absence of molecular features and the presence 
of the Li I 6708 A line, which is typ ically observed in G5- 
K2 giants le.g.. IKumar et aHl2011[) . allows us to assume 
a Teg - 5000 ± 1000 K. 

To put more constraints on the fundamental param¬ 
eters of the cool companion, we compared our spectra 
of MWC 728 with those of Li-rich giants obtained at 
the McDonald Observatory with the same instrumen¬ 
tal setup. The spectra were normalized to the contin¬ 
uum and smoothed with the same box-car filter to get 
rid of a high-frequency noise. A very good match was 
found with the spectrum of the mag netic and Li-rich g i- 
ant HD 232862 (G8 ii, T.ff=5000K. ILebre et all [200^ . 
In the top panels of Figure[6] we compared the spectra 
of MWC 728 and HD 232862 in selected intervals. To fit 
the absorption line profiles of MWC 728, we added a fea¬ 
tureless continuum to the spectrum of HD 232862 in each 
interval. The best agreement was achieved with the cool 
star contribntions of 5%, 7%, 12%, and 15% in the 5410- 
5470 A, 6100-6170 A, 6705-6750 A, and 8320-8390 A 
intervals, respectively. 

The profiles and positions of the metallic absorption 
lines in the spectra of MWC 728 are noticeably vari¬ 
able (Figure[6l bottom panels). Dne to the line weak¬ 
ness that limits the accnracy of individual radial velocity 
measurements, we used the cross-correlation method im¬ 
plemented in the IRAF package rvsao. A wavelength 
region between 6100 A and 6170 A that is well pop¬ 
ulated with the absorption lines was chosen for the 
measurements. The region includes the following lines: 
Ca I (Mult. 3) 6102.72, 6122.22, and 6162.17 A, Fe i 
(Mult. 169) 6136.62 A, Fe i (Mult. 816) 6141.73 A, O i 
(Mult. 10) 6155.99, 6156.78, and 6158.19 A, and Ca i 
(Mult. 20) 6163.76 and 6169.56 A. No Li i line at 6104.00 
A can be clearly separated from the Ca i 6102.72 A line. 
Forty four spectra obtained at CFHT, McDonald, CAN 
SPM, and TCO were used to measure radial velocities. A 
high signal-to-noise spectrum obtained at CFHT on 2004 
December 24 was chosen as a template. The absorption 
lines are well defined in this spectrum which allowed an 
accurate measurement of the averaged radial velocity of 
21.9±1.1 kms“^. All radial velocities measured by cross- 
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Fig. 3. — Left panel: Comparison of the He I 4471 A and the Mg ll 4482 A lines in the spectrum of MWC 728 taken on 2008 December 
15 at McDonald with those of normal stars. The spectrum of BS4967 (B5/6 ill, usini = 115 kms“^) was taken at TCO, and the spectrum 
of BS896 (B6 lll, usinf = 105 kms“^) was taken at the OHP. Intensity is normalized to the nearby continuum, while the spectra were 
moved to the same He I line position for all the stars. Right panel: Tefi relationship with EWs in A of the He 1 4471 A and the Mg ll 
4482 A lines. Circles represent data for normal B—type stars (filled - OHP, open - TCO) and crosses repre s ent OHP data for Be stars. 
Temperature determinations were collected from various papers (e.g.. FF^emat et al.ll2005l : IZorec et al.l[20091 : l^rec &: Hoveil[20131 . The 
horizontal dashed lines show ranges of the EW variations detected in t he spectra of MWC 728, and the vertical dashed lines show most 
likely boundaries of the Tefj for MWC 728 (see explanations in Sect . 13.ill . 


TABLE 3 

Radial velocities of the absorption lines derived by cross-correlation 


JD 

Phase 

RV 

Err. 

JD 

Phase 

RV 

Err. 

JD 

Phase 

RV 

Err. 

JD 

Phase 

RV 

Err. 

3363.730 

0.000 

21.9 

1.1 

4419.892 

0.407 

15.4 

2.1 

6244.761 

0.769 

27.1 

1.1 

6962.795 

0.880 

42.3 

8.1 

3656.801 

0.658 

2.9 

3.5 

4421.886 

0.480 

10.0 

3.4 

6250.694 

0.984 

29.7 

2.2 

6963.764 

0.915 

41.0 

3.7 

3721.654 

0.016 

20.7 

1.9 

4423.866 

0.552 

9.2 

7.8 

6295.634 

0.619 

4.0 

2.1 

6964.758 

0.951 

37.2 

2.6 

3726.711 

0.200 

22.6 

1.9 

4744.841 

0.224 

26.0 

2.3 

6320.678 

0.529 

19.2 

1.7 

6993.739 

0.005 

32.8 

3.7 

4081.875 

0.115 

24.4 

3.1 

4746.855 

0.297 

18.6 

4.3 

6322.655 

0.601 

15.7 

2.7 

6996.699 

0.113 

37.2 

2.9 

4082.658 

0.144 

38.1 

3.2 

4815.747 

0.803 

17.4 

1.8 

6580.981 

0.995 

40.4 

2.2 

7035.618 

0.528 

22.1 

3.6 

4083.884 

0.188 

33.3 

3.7 

5079.925 

0.409 

5.8 

2.6 

6584.811 

0.135 

38.5 

2.8 

7036.733 

0.569 

18.7 

3.8 

4092.612 

0.506 

16.0 

3.5 

5081.934 

0.482 

12.8 

2.7 

6585.786 

0.170 

39.7 

2.3 

7037.609 

0.601 

15.2 

5.8 

4094.576 

0.577 

2.5 

3.1 

5145.750 

0.803 

32.2 

6.3 

6587.831 

0.244 

32.8 

2.3 

7038.613 

0.637 

9.7 

6.3 

4096.750 

0.656 

21.3 

4.7 

5410.634 

0.436 

20.4 

1.6 

6628.678 

0.730 

14.0 

3.5 

7040.720 

0.714 

11.9 

4.3 

4412.847 

0.151 

34.6 

5.0 

5414.580 

0.579 

21.4 

0.8 

6961.798 

0.844 

37.9 

1.9 

7041.717 

0.750 

21.4 

3.8 


The list of Julian dates (2450000-I-) when spectra with signal-to-noise ratios suitable for the absorption line measurements, column 
’’Phase” contains phases calculated for the 27.5-day orbital period, column ”R,V” lists heliocentric radial velocities in kms“^, and 
column ”Err” shows l-cr uncertainties of the radial velocities in kms“^. 


correlation are presented in Table [3l 

The Discrete Fourier Transform (D FT) methocFI 
(iDeemind 1197511 and CLEAN procedure (jR.oberts et all 
Il987ll wCTe applied to search for periodicity in the radial 
velocity data. The resulting periodograms are presented 
in Figure [3 (left panel). The two highest peaks in the 
DFT power spectrum have nearly equal strengths at fre¬ 
quencies that correspond to periods of 29.7 and 27.5 days, 
while the CLEAN procedure picks on the only period of 
Porb= 27.5(1) days. The latter one was adopted as an 
orbital period of the system. The radial velocity curve 
and the E-band photometric data folded on Porb are 
shown in Figure [3 (right panel). Figure[8] shows phased 

^ http://www.univie.ac.at/tops/Period04/ 


time series spectra for the 6100-6170 A wavelength re¬ 
gion where positional variations of the absorption lines 
are clearly seen. The ephemeris of the orbital phase from 
the radial velocity measurements can be described as 

To = HJD 2453356.855 + 27.499(108) x A, 

where Tq corresponds to a moment when the cool com¬ 
panion is located between the hot companion and the 
observer. The radial velocity curve is symmetric about a 
y-velocity of -1-25 kms“^. It can be considered sinusoidal 
within the measurement errors (circular orbit) and has a 
semi-amplitude of ^20 kms“^. 

We note that sometimes the absorption lines are split 
into two components. The splitting develops within a 
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TABLE 4 

Properties oe the Ha line in the spectra oe MWC 728 


MJD 

Iv 

Ir 

EW 

Phase 

Obs. 

MJD 

Iv 

Ir 

EW 

Phase 

Obs. 

3363.730 

6.30 

12.47 

60.2 

0.250 

CFHT 

6661.585 

3.97 

5.33 

26.7 

0.177 

TCO 

3656.801 

3.70 

6.61 

39.5 

0.908 

SPM 

6670.291 

5.49 

5.88 

31.0 

0.493 

Ondrj 

3721.654 

5.85 

6.22 

37.2 

0.266 

McDon 

6670.577 

5.61 

5.97 

31.2 

0.504 

TCO 

3726.711 

5.63 

8.38 

44.9 

0.450 

McDon 

6672.564 

6.22 

7.01 

36.7 

0.576 

TCO 

4081.875 

6.20 

8.06 

50.2 

0.365 

SPM 

6674.564 

5.84 

7.13 

37.7 

0.649 

TCO 

4082.658 

5.21 

8.22 

46.7 

0.394 

SPM 

6676.562 

4.75 

8.08 

36.8 

0.721 

TCO 

4083.884 

5.62 

8.18 

42.8 

0.438 

SPM 

6677.588 

4.52 

8.72 

37.9 

0.759 

TCO 

4092.612 

3.48 

6.82 

39.7 

0.756 

McDon 

6678.553 

3.92 

7.58 

32.1 

0.794 

TCO 

4094.576 

3.86 

7.41 

41.7 

0.827 

McDon 

6680.543 

4.23 

7.58 

33.7 

0.866 

TCO 

4096.750 

4.90 

7.81 

44.4 

0.906 

McDon 

6682.525 

4.86 

7.73 

34.3 

0.938 

TCO 

4412.847 

4.61 

6.65 

40.8 

0.401 

SPM 

6711.328 

5.50 

7.37 

36.7 

0.986 

Ondrj 

4419.892 

3.55 

5.36 

34.2 

0.657 

SPM 

6712.438 

5.60 

7.45 

36.8 

0.026 

Ondrj 

4421.886 

3.80 

5.88 

37.6 

0.730 

SPM 

6712.540 

5.30 

7.86 

37.0 

0.030 

TCO 

4423.866 

4.45 

7.77 

47.9 

0.802 

SPM 

6717.273 

6.06 

6.15 

35.9 

0.202 

Ondrj 

4744.841 

4.71 

6.67 

34.3 

0.474 

SPM 

6861.564 

4.24 

9.22 

43.0 

0.449 

Ondrj 

4746.855 

6.00 

6.49 

33.0 

0.547 

SPM 

6948.702 

4.75 

8.64 

38.7 

0.618 

TCO 

4813.774 

4.75 

8.05 

37.5 

0.981 

McDon 

6949.725 

4.76 

8.41 

37.8 

0.655 

TCO 

4815.747 

5.61 

9.57 

47.5 

0.053 

McDon 

6955.744 

5.09 

7.49 

32.6 

0.874 

TCO 

5079.925 

2.75 

10.37 

44.1 

0.660 

McDon 

6956.777 

5.26 

7.42 

35.5 

0.911 

TCO 

5081.934 

2.62 

8.62 

33.6 

0.733 

McDon 

6957.727 

4.89 

6.69 

33.7 

0.946 

TCO 

5143.723 

6.25 

7.63 

52.4 

0.980 

SPM 

6961.720 

4.98 

11.14 

56.6 

0.091 

TCO 

5145.750 

4.71 

8.82 

50.2 

0.053 

SPM 

6961.798 

5.28 

12.05 

61.9 

0.094 

SPM 

5410.634 

5.54 

6.35 

34.2 

0.686 

CFHT 

6962.795 

4.30 

11.67 

57.1 

0.130 

SPM 

5414.580 

5.64 

6.35 

34.1 

0.829 

CFHT 

6963.764 

3.67 

11.16 

50.3 

0.166 

SPM 

5925.657 

4.90 

6.60 

29.5 

0.415 

TCO 

6964.654 

3.23 

11.82 

47.3 

0.198 

TCO 

6154.576 

2.84 

7.05 

29.2 

0.739 

Ondrj 

6964.758 

3.42 

12.03 

49.6 

0.202 

SPM 

6158.531 

2.57 

5.91 

25.7 

0.883 

Ondrj 

6969.693 

2.16 

8.77 

31.2 

0.381 

TCO 

6187.475 

4.40 

5.56 

27.8 

0.936 

Ondrj 

6971.685 

3.39 

7.71 

29.4 

0.454 

TCO 

6244.761 

5.16 

7.89 

44.1 

0.019 

SPM 

6972.701 

3.80 

7.18 

28.8 

0.491 

TCO 

6250.694 

4.20 

5.66 

30.9 

0.235 

SPM 

6973.731 

4.07 

6.60 

28.3 

0.528 

TCO 

6290.675 

4.02 

7.42 

23.3 

0.689 

McDon 

6976.659 

4.24 

6.22 

30.4 

0.635 

TCO 

6295.634 

6.13 

5.71 

38.3 

0.869 

McDon 

6980.685 

4.47 

5.72 

35.0 

0.781 

TCO 

6320.678 

3.45 

5.58 

22.1 

0.780 

SPM 

6982.656 

4.32 

6.43 

33.9 

0.853 

TCO 

6322.655 

2.59 

3.57 

14.7 

0.852 

SPM 

6992.652 

4.74 

5.85 

29.7 

0.216 

TCO 

6580.981 

4.83 

4.54 

33.7 

0.246 

SPM 

6993.718 

4.62 

6.61 

36.5 

0.255 

SPM 

6584.811 

4.84 

8.76 

49.0 

0.385 

SPM 

6996.678 

5.23 

7.85 

44.5 

0.362 

SPM 

6585.786 

4.45 

9.19 

49.3 

0.420 

SPM 

7002.668 

2.89 

6.42 

24.1 

0.580 

TCO 

6587.831 

1.93 

10.30 

38.4 

0.495 

SPM 

7003.596 

2.64 

6.56 

24.0 

0.614 

TCO 

6589.851 

1.97 

8.92 

28.2 

0.568 

SPM 

7004.645 

2.95 

6.30 

23.8 

0.652 

TCO 

6591.920 

2.80 

7.40 

25.1 

0.643 

SPM 

7006.590 

3.50 

5.90 

22.5 

0.723 

TCO 

6592.607 

2.49 

6.63 

21.4 

0.668 

Ondrj 

7032.592 

3.47 

9.88 

42.6 

0.669 

TCO 

6592.956 

2.97 

7.89 

24.8 

0.681 

SPM 

7033.570 

2.62 

10.28 

38.8 

0.704 

TCO 

6600.758 

4.88 

7.74 

34.7 

0.965 

TCO 

7035.618 

2.04 

9.44 

32.9 

0.779 

SPM 

6605.735 

4.26 

7.95 

38.3 

0.146 

TCO 

7036.733 

2.56 

9.02 

31.6 

0.819 

SPM 

6606.506 

3.74 

7.15 

32.0 

0.174 

Ondrj 

7037.609 

2.44 

8.39 

29.2 

0.851 

SPM 

6607.710 

4.18 

7.49 

33.7 

0.218 

TCO 

7038.613 

2.57 

7.63 

26.8 

0.887 

SPM 

6611.709 

4.88 

6.49 

30.4 

0.363 

TCO 

7039.559 

3.28 

7.29 

25.6 

0.922 

TCO 

6625.712 

3.59 

8.03 

34.8 

0.872 

TCO 

7039.748 

3.48 

7.28 

26.4 

0.929 

SPM 

6626.663 

4.20 

7.21 

35.0 

0.907 

TCO 

7040.720 

3.43 

6.35 

26.5 

0.964 

SPM 

6627.690 

4.22 

6.83 

31.3 

0.944 

TCO 

7041.586 

3.14 

6.65 

27.4 

0.996 

TCO 

6628.678 

4.33 

6.52 

30.1 

0.980 

TCO 

7041.717 

3.25 

6.51 

27.4 

0.000 

SPM 

6637.655 

3.46 

6.61 

31.9 

0.307 

TCO 

7042.595 

3.58 

6.45 

27.4 

0.032 

TCO 

6639.658 

3.67 

7.62 

36.6 

0.379 

TCO 

7044.508 

4.24 

5.89 

25.9 

0.102 

TCO 

6642.650 

3.90 

6.84 

32.3 

0.488 

TCO 

7045.529 

4.57 

5.90 

26.5 

0.139 

TCO 

6643.691 

3.81 

7.18 

33.4 

0.526 

TCO 

7059.535 

5.96 

5.71 

36.5 

0.648 

TCO 

6645.682 

4.06 

5.47 

30.4 

0.598 

TCO 

7061.549 

4.98 

6.52 

42.3 

0.722 

TCO 

6653.641 

7.27 

7.69 

52.0 

0.888 

TCO 

7062.544 

5.68 

8.14 

47.2 

0.758 

TCO 

6654.627 

5.83 

8.32 

50.4 

0.924 

TCO 

7065.554 

6.32 

8.44 

50.5 

0.867 

TCO 


The columns lists Modified Julian dates (JD—245000), the blue (Iv) and red (Ij^) peak strengths in the continuum units, the Ha 
line EW in A, orbital phase according to the ephemeris in the text, and the observing site ID (see Table[T]l, respectively. 
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Wavelength[A] 


Fig. 4. — Comparison of the blue part of the spectrum of 
MWC 728 with those of BS896 (B6 iii) and Vega (AO v). The 
spectra of Vega and MWC 728 were obtained at OAN SPM, while 
the spectrum of BS 896 was taken at OHP (4000 A is its shortest 
wavelength). 


I/Ic 



5860 






5880 5890 

Wavelength [A] 


5900 


Fig. 5. — Variations of the He i 5876 A line. The spectra were 
obtained (from bottom to top) at: CFHT in 2004, the McDon¬ 
ald Observatory in 2008 and 2012, and OAN SPM in 2014. The 
wavelength scale is heliocentric. The vertical dashed line shows the 
rest wavelength of the line at the systemic radial velocity of +25 
kms“^. Individual spectra are shifted by 0.15Ic with respect to 
each other. Numerous telluric lines are clearly visible around the 
interstellar Na i D—lines (5889 and 5895 A), especially in the two 
bottom spectra. The sodium lines have stable positions and are 
shown here to contrast the He line variations and demonstrate the 
wavelength calibration accuracy. 


few days and seems to occur at no preferential orbital 
phase. The components appear on both sides of the ini¬ 
tial unsplit line position and are separated by up to 80 
kms“^. Such an episode observed on 2006 December 22~ 
26 is shown in Figure[S] (bottom panels). Splitting occurs 
simultaneously in all absorption lines which adds to the 
scatter of the radial velocities determined by cross corre¬ 
lation. This phenomenon might be due to appearance of 


emission components near the line centers, but its actual 
mechanism is unclear. 

The third group of absorption lines has interstellar 
origin. MWC 728 is located at a galactic latitude of 
—20°, nearly 3° south of the association PerOB2. 
The distance of Per OB 2 from the Sun is 0.3 kpc (e.g., 
iMel’nik fc PambT^ 1200911 . and the inter stellar extinc tion 
in its line of sight is Ay = 0.95 mag (|Cernisl[l993[l . A 
photome t ric stu dy of a ~ 1° area around MWC 728 by 
IRvdgrenI (|1971ll showed that the extinction in this direc¬ 
tion is lower and reaches 1 mag at a distance of ~1 kpc. 
The interstellar absorption lines in the MWC 728 spec¬ 
trum are represented by diffuse interstellar bands (DIBs), 
such as 4430, 5780, 5797, 6613 A, and the Na i D-lines 
at 5889 and 5895 A. The sodium lines are narrow and 
saturated at our highest resolution spectra (R ^60000). 
Their measured average radial velocity is 11±1 kms“^. 
The DIBs are weak (e.g., EW 5780 = 0.13 A). Based 
on the relatio nship of the EW 5780 and E(B — V) from 
iHerbid (1199311 . the interstellar reddening for MWC 728 
is E{B — V) ~ 0.3 mag. This estimate implies a visual 
interstellar extinction of Ay ^ 1.0 mag, and a distance 
to the object of at least 1 kpc. Both the distance and 
reddening are discussed in more detail in Sect.[4l 

3.2. Emission lines 

The emission-line spectrum is represented by the 
Balmer lines, the [O i] 6300 A and 6364 A lines (Eig- 
ure ED, and weak emission components of the He i 5876 A 
line mentioned above (Figure E]). The forbidden oxygen 
lines have double-peaked profiles with nearly the same, 
slightly varying intensity peaks and an average peak sep¬ 
aration of 40 kms“^. The average radial velocity of these 
lines is +27 ± 5 kms“^ which virtually coincides with 
the 7 -velocity of the absorption lines. The equivalent 
widths of the [O i] lines are EWeaoo = —0.9 ± 0.1 A and 
EW6364 = -0.34 + 0.05A. 

The Ha and H/3 lines are relatively strong and show 
double-peaked emission profiles in most of our spectra 
(Figure [9] and [To]). Measured parameters of the Ha line 
are presented in Table|4| and plotted in Figuredl] (for 
the seasons of 2013-2014). The Ha line profile is highly 
variable with changes occurring on a time scale from days 
(Figure fTUl) to months iFigurc fTTl) . The EW of the Ha 
line emission component varies from —15 A to —61 A. 

The blue peak of the Ha line is usually weaker than 
the red one with an average peak intensity ratio of V/R 
= 0.6+0.2. This is different from the Balmer line profiles 
of classical Be stars, in which the peaks typically show 
a nearly equal strength (except for cyclic variations of 
the peak ratio, V/R variations, that are due to density 
perturbations rotating in the disk). This effect may be 
explained by emission from the stellar wind of the hot 
companion that expands beyond its disk (see Eigurc[T3|). 

Sometimes the peak intensity ratio reaches V/R ~ 1 
(Eig|T2|). The Ha line EW is typically —35 + 3 A at these 
moments. However, on one occasion (2013 December 26) 
the line was stronger with an EW of —52 A. We interpret 
the profiles with V/R « 1 as dominated by the circum- 
primary disk emission with almost no contribution from 
the primary’s stellar wind. The EW changes of the sym¬ 
metric Ha profiles may be due to density variations in 
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Fig. 6.— Top panels: Absorption lines in the McDonald spectra of MWC 728 and HD 232862 (G8 II, see text) for three spectral regions. 
Lines of neutral metals are marked. Bottom panels: Examples of the short-term variability and splitting of some absorption lines in 
different regions of the MWC 728 spectra are presented. Vertical lines show laboratory positions of the Ca I, Li I, and Ti I. 




Fig. 7.— Left panels. Power spectra calculated by the DFT (bottom) and CLEAN (top) algorithms for the absorption lines radial 
velocity derived by cross-correlation (see text). Right panel. Radial velocity (top) and the V-band photometry (bottom) folded on the 
proposed orbital period Porb= 27.5 days. The V-band brightness is given with respect to that of the comparison star, HD 281193. 


the circum-primary disk. 

We do not hnd noticeable variations of the Ha line 
profile related to the periodicity of the absorption lines. 
However some imprints of the orbital motion can be seen 
in the violet peak in Figure d^ where phased time series 
spectra around Ha folded with the orbital period are 
shown. The current amount of data is not sufficient to 
verify this suggestion. A few months of continuous spec¬ 
troscopic observations are needed to explore the validity 
of this statement. 

Another feature of the Ha profile is a weak emission 
peak at a nearly zero heliocentric radial velocity. It is 


clearly visible in most of our high-resolution spectra but 
occasionally overlaps with one of the main peaks (Figs.[9l 
[ini and [HD. The source of this emission is not clear yet. 

To summarize our current understanding of the forma¬ 
tion and variations of the Ha line we roughly reproduced 
it by a simple toy model (Figure [13]). It includes an at¬ 
mospheric absorption from the hot star, a double-peaked 
emission from the circum-primary disk, an emission from 
the shell/wind with a P Gyg-type profile, and a weak 
emission peak at the systemic velocity. All the com¬ 
ponents were represented by simple Gaussian profiles. 
The resulting profile was calculated following equation 
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Fig. 8.— Phased time series spectra in the 6100 - 6170 A wavelength region folded three times (Y-axes) on the orbital period of the 
system Poj.b= 27.5 days. The intensity scale in the continuum units is shown in the right part of the graph. 



Wavelength fAl Wavelength fAl 

Fig. 9. — Typical profiles of the H/3, Ho; and [O i] emission 
lines in the spectrum of MWC 728. The spectrum was taken on 
10/30/2014 at OAN SPM. The wavelength scale is heliocentric, the 
intensity is normalized to the nearby continuum. 


12 from iSchoenberd ()1985h as /a = S / x I'pcyg~^^Wcyg^ 
where E / is a sum of all normalized fluxes of the under¬ 
lying components and I'pcygi ^'pcyg absorption and 
emission components from the P Cyg-type profile of the 
shell/wind. Such a toy model can qualitatively repro¬ 
duce the observed profiles and their variability. However 
application of a radiation transfer code that is capable of 
treating variable circum-primary disk and wind compo¬ 
nents is required for quantitative modeling. The latter is 
beyond the scope of this paper. 

3.3. Brightness variations 

Photometric observations of MWC 728 are sparse. It 
was observed in the co urse of the Northern Stars Vari¬ 
ability Survey (NSVS, iWozniak et al.l 12004( 1 in 1999- 
2000. The object’s brightness was detected to be 
V=10.17±0.05 mag with no correction to the standard 
Johnson system. No regular variations or evidence for 
eclipses were found in these data. 

Our photometric data obtained at TShAO show a scat¬ 
ter of ^0.1 mag (see Figure[7]), similar to that of the 
NSVS data. However on two nights (2014 November 2 
/ JD2456964 and 2015 February 4-5 / JD2457058) we 
managed to take photometry contemporaneously with 
times when the emission-line spectrum got stronger. On 
both these occasions the EWhu increased by a factor of 


^1.5 from ^ —40 A to ~ —60 A, but the object faded by 
~ 0.1—0.25 mag in the V-band from the quiescence level. 
These variations of both the spectral line and brightness 
can be explained by sporadic matter outbursts from the 
surface of the hot star. Part of this material obscures the 
stellar surface and leads to a fading, while line emission 
becomes stronger due to increasing amount of circumstel- 
lar gas. A similar phenomenon has been observed during 
the built-up of the circ um-primary disk around in the S 
Scorpii binary system ((Miroshnichenko et al.l[200^ . 

4. DISCUSSION 


4.1. The binary system parameters 

Adopting an orbital period of 27.5 days and a semi¬ 
amplitude of the radial velocity curve of K 2 = 20 kms“^, 
we derive a mass function of 


PorbATl 

27rG 


Ml sin^ i 

{l + qf 


^ 2.3 X IO-^Mq, 


where Mi is the primary (hot) companion mass, q = 
M 2 /M 1 is the secondary to primary mass ratio, and i is 
the orbital inclination angle. 

The small mass function implies a low inclination an¬ 
gle to be reconciled with a meaningful hot star mass and 
the absence of the photometric evidence for eclipses. No 
signs of the orbital motion in the Balmer lines formed 
in the circum-primary disk require a g < 1 (see Fig- 
urc ITSl) . Additionally, the detection of the radial velocity 
variations of the absorption lines and the double-peaked 
Balmer emission-line’s profiles are inconsistent with the 
system viewed exactly pole-on. 

As mentioned in Sect. 13.11 the He i 4471 A and Mg ii 
4482 A line profiles suggest a wsinz = 110 kms“^ for 
the hot star (Figure [31 left panel). Single main sequence 
stars with spectral types B5±l have masses of 3.5Mq < 
Ml < 4 . 5 M 0 , radii of ^3.3 R fA, and breakup rota tional 
velocities of ^^420 kms“^ fe.g. JFremat et al.ll2005f l. The 
presence of a circum-primary disk suggests that the hot 
star rotates near the breakup velocity, and an upper limit 
on the orbital inclination angle is ~15° if the circum- 
primary disk and the orbital plane are coplanar. A more 
massive primary would require an z < 10° (Figure [15]) • 
Assuming g = 1 as an upper limit for the mass ratio leads 
to the companions’ separation of < 74R0 and the size of 
the secondary’s Roche lobe of <26Ro. 

There is a potential problem in explanation of the 
Balmer line profiles whose relatively deep central depres¬ 
sions (see, e.g., Figurc[4| and [9]) may not be consistent 
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Fig. 10.— Short-term spectral variations of the Ho (left panels) and H/3 (right panels) lines in the spectrum of MWC728 in October 
2013 and November 2014. The radial velocity scale is heliocentric, the intensity is normalized to the nearby continuum. 
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JD-2450000 JD-2450000 


Fig. 11.— Variations of the H« line in the spectrum of MWC 728 detected during the 2013 (left panel) and 2014 (right panel) observing 
seasons. From bottom to top are shown: equivalent width (EW) in A, blue-to-red peak intensity ratio (V/R), heliocentric radial velocity 
of the red (Vi.ed) and blue (Vniue) emission peak in kms“^. 


with the derived low orbital inclination angle. This ap¬ 
parent contradiction may be resolved by assuming pre¬ 
cession of the circum-primary disk. Such an explanation 
has already been proposed to account for emission -line 
profile variations in some Be stars ('e.g.. lHiratal[200l . In 
the case of MWC 728, the precession period should be 
large as we did not detect any noticeable changes in the 
emission-line profiles for over 10 years. 


4.2. Relative flux contributions 

In order to place more tight constraints on the funda¬ 
mental parameters of the companions and the distance 
to the system we need to estimate relative contributions 
of all the system constituents, i.e. both stars and the 
circumstellar matter. The cool companion contribution 
cannot be determined precisely due to veiling by the 
other flux components that prevents us from learning 
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Fig. 12.— The profiles of the Ho lines when V/R ratio is close to 
one. The spectra are normalized to continuum and shifted along 
the Y-axis for the best presentation. The radial velocity scale 
is heliocentric.The vertical dashed line shows the systemic radial 
velocity of +25 kms~^. 



Fig. 13.— A toy model of the Ha profile of MWC 728 observed 
on 10/30/2014 at OAN SPM (see description in the text). The 
intensity and heliocentric radial velocity are shown in the units as 
in Figure[To] 

its true rotation rate and therefore the strength of its 
spectral lines. Taking this uncertainty into account and 
using our results for the spectral line contributions (see 
Sect. l3.1ll . we adopt a l/-band contribution of 10% from 
the cool companion to the total syste m flux. 

Using stellar atmosphere models of iKuruci (Il9^ , we 
found that a 5000 K star continuum contributions to 
that of a 14000 K star are nearly the same as the spec¬ 
tral line contributions in all the wavelength intervals we 
mentioned in Sect. EH). This implies that the disk con¬ 
tribution to the total flux has nearly the same wave¬ 
length dependence as that of the hot star at least be¬ 
tween A ^^5390 A and A 8800 A. A very similar result 
has been derived for optical e mission from disks of Be 
stars le.g.. iCarciofi et al1l2006ll . 

Since most of the disk surface seems to be facing the 
observer, it should result in a non-zero disk contribution 


to the system radiation. The J-band excess is mostly 
due to the free-free and bound-free disk radiation, be¬ 
cause circumstellar dust typically has temperatures be¬ 
low 1500 K and dominates the spectral energy distribu¬ 
tion at longer wavelengths (see FigureE]). At the same 
time, we can estimate upper limits on the J-band excess 
and interstellar extinction by assuming no circumstellar 
contribution to the continuum in the optical range. Us¬ 
ing our average photometric data for the dates when no 
strong spectral variations were observed (U = 9.77 mag, 
B — V = 0.27 mag, see Table[2]) and taking into account 
the contribution from the cool companion, we get the fol¬ 
lowing upper limits: E{B — V) = 0.39 mag. Ay = 1.22 
mag, and the J-band excess of ^0.3 mag. 

Based on these results, we adopt a lower interstel¬ 
lar extinction of Ay 1.0 mag, which is consistent 
with that derived from the diffuse interstellar band EW 
(see Sect. 13.II) and implies a circumstellar color-excess of 
E{B — V) ^0.1 mag. This estimate does not contradict 
our observations of the He i 4471 A and Mg ii 4482 A 
lines and suggests a disk contribution of <20 % to the 
continuum at these wavelengths. 

Finally, we estimate the expected contributions of the 
hot companion, the disk, and the cool companion in the 
optical range (in the U-band) to be roughly 60%, 30%, 
and 10%, respectively. This is in agreement with a gen¬ 
eral conclusion that the disk contribution to the total op¬ 
tical flux in Be stars may hard l y exce ed 50% of the hot 
companion flux iHaubois et al.l (120121) . The wavelength 
dependence of the disk radiation in the IR region was 
calculated by extrapolation of the difference between the 
total flux and that of both stars (see Fig[2|). 

The fractional contributions define a lower limit the 
secondary’s radius R-? = 7.7Ri^ (correspondi ng to a lu¬ 
minosity type III, iStraizvs fc Kurilienelll98ll ). if we as¬ 
sume a typical main-sequence mass of Mi = A.OMq and 
a radius of Ri = 3.3i?0 for the primary. The compan¬ 
ions would have these radii at a distance of 1 kpc, which 
is consistent with the distance dependence of the inter¬ 
stellar extinction (see Sect l3.ip and can be considered a 
lower limit of the system distance. An upper limit is 
^3.5 kpc at which the secondary fills its Roche lobe. 

Fundamental parameters of the companions for both 
distances are listed in Table [5) The upper distance limit 
is hardly consistent with the orbital parameters, because 
it corresponds to a low critical rotational velocity of the 
hot star (^320 kms“^). The latter requires i > 20° and 
a cool companion that is more massive than the hot star 
(see FigurefTSl). Therefore, it is unlikely that the cool 
companion fills its Roche lobe. This allows us to assume 
that the distance to the system is close to its lower limit, 
and the companions mass ratio is q ^0.5. 

4.3. Possible system model 

The strong IR excess at A > 2^m (see Figured]) im¬ 
plies the presence of dust in the system. It should 
exist in the circum-binary area, as the stars are too 
close together to allow its existence around either star. 
iMiroshnichenkol (|2007f ) proposed that the dust forms out 
of the material lost from the system in the course of 
mass-transfer between the stars. Evolution of interme- 
di ate-mass binary system s has b een studied theoretically 
bv ivan Rensbergen et al.l (|2008D . They calculated a grid 
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Fig. 14.— Phased time series spectra around the Ha line folded three times (Y-axis) on the orbital period of the system Poj.i 3 = 27.5 day. 
The color coded intensity scale in the continuum units is shown on the right. 



Fig. 15.— M 2 vs Ml relationship for different system inclina¬ 
tions. The dashed lines correspond to mass ratios q = 0.5 and q 
= 1. The shaded region marks the proposed mass range of the hot 
star based on its spectral type. The part of this region below the 
q = 1 line marks the possible range of masses of the secondary. 


TABLE 5 

Fundamental parameters of the stellar companions of the 
MWC 728 BINARY system at different distances 


D, kpc 

logLi/Le 

Ri/Rq 

logL2/T0 

R 2 /RQ 

1.0 

2.8 

3.3 

1.6 

7.7 

3.5 

3.7 

15.0 

2.6 

26.0 



log Teff 


Fig. 16.— Evolutionary tracks of a 5 Mq (solid line) and a 2 
Mq (dotted line) compone nts of a close binary system with non¬ 
conservative mass loss from ivan Rensbergen et al.l II2008I1 . The on¬ 
set of mass transfer is marked with an open circle on the track of 
the 5 M0 star, while the 2 Mq star is still near the starting point 
of its evolution. The tracks are shown until the end of the mass 
transfer stage. The dashed lines show the zero-ag e and terminal- 
age m ain sequence for rotating single star models (lEkstrom et al.l 
1201211 . The filled circles show the fundamental parameters of both 
companions for the lower limit on the system distance, and the 
filled squares show them for the higher distance limit (see Tablel^ . 


Column 1 lists the distance, columns 2 and 3 show the lumi¬ 
nosity and radius of the hot star, columns 4 and 5 show the 
luminosity and radius of the cool star. We assumed Tgfi = 
14000 K and BCrz = —1.1 mag for the hot star, Teff = 5000 
K and BCi/ = —0.3 mag for the cool star. The b olometric 
corrections are taken from IMiroshnichenkol 119971 1. 

of models for conservative evolution as well as for a weak 
and strong tidal interaction between the stars. The 
grid includes detailed calculations of physical and or¬ 
bital parameters from formation through the entire mass- 
transfer phase and beyond. Non-conservative (liberal) 
evolution, when some mass is lost from the system, was 
only considered for the more massive companion mass of 
>5 Mq. 

Browsing through the grid, we found a model which 
qualitatively explains our findings for MWC 728. It has 
initial masses of 5 Mq and 2 Mq and an initial orbital 
period of 8 days. The mass-transfer starts in 1.1210® 


years since the evolution began when the more massive 
star fills its Roche lobe. It lasts for ~ 5.5 10® years when 
the donor reaches a mass of 0.9 Mq, the gainer has 4.2 
Mq , and nearly 2 Mq is lost from the system. Therefore 
there is enough material to produce dust in the circum- 
binary area. The evolutionary tracks for the stars in this 
model are shown in Figure flHl 
Although the orbital period evolution in the model 
does not match that of MWC 728 (it reaches 27 days 
in the model, but only when the companions have al¬ 
most the same temperature and luminosity), it is qualita¬ 
tively consistent with the nature of this system. In parti¬ 
cular, the companions’ luminosity and their ratio along 
with their temperatures roughly match our results. Re¬ 
fining the system parameters with further observations 
and calculating a finer grid of binary models with non¬ 
conservative evolution would provide a more complete 
explanation of the system properties. 
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5. CONCLUSIONS 

Our study confir m an earlier suggestion by 
iMiroshnichenko et alJ (|2007t) that MWC 728 is a 
binary system. We found the orbit to be circular (within 
the measurement errors) with a period of Porb = 27.5 
days. The system consists of a B5 primary that is '^6 
times brighter (in the U-band) than a G8 ill secondary 
and is most likely located at a distance of kpc. 

The presented results further support the idea that 
FS CMa objects are a group of moderately evolved binary 
systems and not proto-planetary nebulae or young stars. 
Although nearly a dozen FS CMa objects shows signs of 
binarity mo stly through the presence of spectral lines of 
a cool star (jMiroshnichenkcill2007t IMiroshnichenko et al.l 
120071 l2011h . MWC 728 is only the second group object 
with detected orbit al motion. The first on e was Cl Cam 
(Porb = 19.41 davs. lBarsukova et al.ll200^ . and it seems 
to have a degenerate secondary companion. In all recog¬ 
nized FS CMa type binaries but one (MWC 623. IZickgrall 
l2001tl spectral lines of the secondary companion are very 
weak and require both high-resolution and high signal-to- 
noise spectroscopy to detect them. Additionally, irregu¬ 
lar variations of the emission-line spectrum contribute 
to veiling of the hot star orbital motion. This phe¬ 
nomenon occurs on a time scale of a few days and may 
represent fast variations of the hot star wind specific to 
this evolutionary stage. These findings call for frequent 
and long-term coordinated spectroscopic and photomet¬ 
ric campaigns to reve al the objects’ natu re. _ 

A recent result by Ide la Fuente et al] (|2015f) on two 
newly found FS CMa group candidates in young clusters 
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brought up an idea that some FS CMa objects could re¬ 
sult from a merger in a binary system. This is definitely 
a possibility, although failure to detect of signs of the sec¬ 
ondary companion does not unequivocally imply its ab¬ 
sence. Models of non-conservative binary evolution (see 
Sect.S]) shows that the cooler secondary may be much 
fainter than the hot primary for a long time. In any 
case, the binary model remains the major mechanism to 
explain objects of this group. 

We will continue coordinated spectroscopic and pho¬ 
tometric observations of MWC 728 in order to further 
study properties of occasional increases of the emission¬ 
line spectrum. They seem to be associated with the ob¬ 
ject’s brightness decreases and should help to better con¬ 
strain the nature of the binary system. 
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